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We analytically describe the architecture of randomly damaged uncorrelated networks as a set of 
successively enclosed substructures — fc-cores. The fc-core is the largest subgraph where vertices 
have at least fc interconnections. We find the structure of fc-cores, their sizes, and their birth points 
— the bootstrap percolation thresholds. We show that in networks with a finite mean number Z2 
of the second-nearest neighbors, the emergence of a fc-core is a hybrid phase transition. In contrast, 
if Z2 diverges, the networks contain an infinite sequence of fc-cores which are ultra-robust against 
random damage. 
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Introduction. — Extracting and indexing highly inter- 
connected parts of complex networks — communities, 
cliques, cores, etc. — as well as finding relations between 
these substructures is an issue of topical interest in net- 
work research, see, e.g., Refs. 0,0. This decomposition 
helps one to describe the complex topologies of real- world 
networks. In this respect, the notion of fc-core is of fun- 
damental importance [3|, |^ . The fc-core may be obtained 
in the following way. Remove from a graph all vertices of 
degree less than fc. Some of the rest vertices may remain 
with less than fc edges. Then remove these vertices, and 
so on until no further removal is possible. The result, 
if it exists, is the fc-core. Thus, a network is organized 
as a set of successively enclosed fc-cores, similarly to a 
Russian nesting doll. 

The fc-core decomposition was recently applied to a 
number of real- world networks (the Internet, the WWW, 
cellular networks, etc.) |a,|g,L| and was turned out to be 
an important tool for visualization of complex networks 
and interpretation of cooperative processes in them. Rich 
fc-core architectures of real networks was revealed. Fur- 
thermore, a fc-core related Jellyfish model Q is one of 
the popular models of the Autonomous System graph 
of the Internet. The notion of the fc-core is a natu- 
ral generalization of the giant connected component in 
the ordinary percolation [53, Hoi Hl| (for another possible 
generalization, see clique percolation in Ref. |l2j1. Im- 
pressively, the giant connected component of an infinite 
network with a heavy-tailed degree distribution is robust 
against random damage of the net. The fc-core perco- 
lation implies the emergence of a giant fc-core below a 
threshold concentration of vertices or edges removed at 
random. In physics, the fc-core percolation (bootstrap 
percolation) on the Bethe lattice was introduced in Ref. 
[llj for describing some magnetic materials. Note that 
the fc > 3-core percolation is an unusual, hybrid phase 
transition with a jump of the order parameter as at a first 
order phase transition but also with strong critical fluc- 
tuations as at a continuous phase transition [HEl- The 
fc-core decomposition of a random graph was formulated 
as a mathematical problem in Refs. 0,0]. This attracted 



much attention of mathematicians 0, ^| , but actually 
only the criteria of emergence of fc-cores in basic random 
networks were found. 

In this Letter we derive exact equations describing the 
fc-core organization of a randomly damaged uncorrelated 
network with an arbitrary degree distribution. This al- 
lows us to obtain the sizes and other structural charac- 
teristics of fc-cores in a variety of damaged and undam- 
aged random networks and find the nature of the fc-core 
percolation in complex networks. We apply our general 
results to the classical random graphs and to scale-free 
networks, in particular, to empirical router-level Internet 
maps. We find that not only the giant connected compo- 
nents in infinite networks with slowly decreasing degree 
distributions are resilient against random damage, as was 
known, but their entire fc-core architectures are robust. 

Basic equations. — We consider an uncorrelated 
network — a maximally random graph with a given 
degree distribution P(q) — the so-called configuration 
model. We assume that a fraction Q = 1 — p of the 
vertices in this network are removed at random. The 
fc-core extracting procedure results in the structure of 
the network with a fc-core depicted in Fig. 

Taking into account the tree-like structure of the in- 
finite sparse configuration model shows that the fc-core 
coincides with the infinite (fc— l)-ary subtree [13 • (The 
m-ary tree is a tree, where all vertices have branching at 




FIG. 1: The structure of a network with a fc-core. The fc- 
core is the internal circle. Vertices of degrees smaller than fc 
form the light-grey area. The dark-grey regions are numerous 
finite clusters with those vertices of degrees q > fc, which do 
not belong to the fc-core. These clusters are either connected 
to the fc-core by less than fc edges or isolated from it. 
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least m.) Let R be the probability that a given end of an 
edge of a network is not the root of an infinite (k— l)-ary 
subtree. Then a vertex belongs to the fc-core if at least fc 
its neighbors are roots of infinite (fc— l)-ary subtrees. So 
the probability that a vertex is in the fc-core is 

M(k) = P J2P(q) ]T Cl&- n {l R) n , (1) 



n— k 



where C™ = m!/(m — n)!n!. Note that for the ordinary 
percolation we must set k = 1 in this equation. 

An end of an edge is not a root of an infinite (fc— l)-ary 
subtree if at most k— 2 its children branches are roots of 
infinite (fc— l)-ary subtrees. This leads to the following 
equation for R: 



fc-2 



n=0 



(t+ l)P(t+l) 
21 



G*iJ'- n (l-ii) n 



(2) 



Let us explain this equation, (i) The first term, 1—p = Q, 
is the probability that the end of the edge is unoccupied, 
(ii) C l n R l ~ n {l — R) n is the probability that if a given 
end of the edge has i children (i.e., other edges than the 
starting edge), then exactly n of them are roots of infi- 
nite (k— l)-ary subtrees. (i+l)P(i+l)/z\ is the proba- 
bility that a randomly chosen edge leads to a vertex with 
branching i. zi = J2 q l^il) i s the mean number of the 
nearest neighbors of a vertex in the graph. Thus, in the 
square brackets, we present the probability that a given 
end of the edge has exactly n edges, which are roots of 
infinite (fc — l)-ary subtrees, (iii) Finally, we take into 
account that n must be at most k — 2. 

The sum X)n=o m © may be rewritten as: 



fe-2 



(l-R) n d n 
n\ dR> 



-Gi(R), 



(3) 



where G\(x) 



i 1 Y uj P(q)qx' 1 - 1 = z^ 1 dG (x)/dx, and 
)x q |ls|. Then Eq. © takes the form: 



R= l-p + p$ k (R). (4) 

In the case p = 1, Eq. was recently obtained in (lif . 
If Eq. has only the trivial solution R = 1, there is 
no giant fc-core. The emergence of a nontrivial solution 
corresponds to the birth of the giant fc-core. It is the 
lowest nontrivial solution R< 1 that describes the fc-core. 
Let us define a function 



f k (R) = [1 - * fc (iJ)]/(l - R). 



(•5) 



This function is positive in the range R € [0, 1) and, 
in networks with a finite mean number of the second 
neighbors of a vertex, Z2 = J2 q l{l — l)P(l)> ^ tends to 
zero in the limit R — > 1 as fk(R) oc (1 — R) k ~ 2 . In terms 
of the function fk(R), Eq. Q is especially simple: 



pfk(R) = 1. 



(6) 



Depending on with increasing R, fk(R) either (i) 

monotonously decreases from /fe(0) < 1 to /fe(l) = 0, 
or (ii) at first increases, then approaches a maximum 
at i? max G (0,1), and finally tends to zero at R — > 1. 
Therefore Eq. © has a non-trivial solution R < 1 if 



p max fk(R) > 1. 
Re[o,i) 



(7) 



This is the criterion for the emergence of the giant k- 
core in a randomly damaged uncorrelated network. The 
equality in Eq. J7J) takes place at a critical concentration 
p c (k) when the line y(R) = l/p c (k) touches the max- 
imum of fk(R)- Therefore the threshold of the fc-core 
percolation is determined by two equations: 



Pc(k) - V/fcCRn 



— /fc(-Rmax)- 



(8) 



-Rmax is the value of the order parameter at the birth 
point of the fc-core. At p < p c (k) there is only the trivial 
solution R = 1. 

At k = 2, Eq. describes the ordinary percolation 
in a random uncorrelated graph [ToL Hll |. In this case, in 
infinite networks we have i? m ax — * 1, and the criterion 
(0) is reduced to the standard condition for existence of 
the giant connected component: pG'^l) = pz2/z\ 1. 

Let us find R near the k > 3-core percolation transition 
in a network with a finite Zi. We examine Eq. for 
R = i? max + t and p = p c (k) + e with e, |r| <C 1. Note 
that at k ^ 3, ^(-R) is an analytical function in the range 
R £ [0, 1). It means that the expansion of $>f-(R+r) over 
r contains no singular term at R € [0,1). Substituting 
this expansion into Eq. Q, in the leading order, we find 



Rn 



-Rcxlp-Pcik)] 1 / 2 , 



(9) 



i.e., the combination of a jump and the square root crit- 
ical singularity. The origin of this singularity is an in- 
triguing problem of the hybrid phase transition. 

The structure of the fc-core is essentially determined 
by its degree distribution which we find to be 

Pk(q) = E P{d')CiR q '- q {l - R) q - (10) 



The mean degree of the fc-core vertices is z\{k) = 
J2 q >k Pk(o)Q- The fc-core of a given graph contains the 
fc+ 1-core as a subgraph. Vertices which belong to the fc- 
core, but do not belong to the k + 1-core, form the fc-shell 
of the relative size S(k) = M(k) - M(k + 1). 

We apply our general results to two basic networks. 

Erdos-Renyi (ER) graphs. — These random 
graphs have the Poisson degree distribution 
P{q) — Ziexp(—Zi)/q\, where Z\ is the mean de- 
gree. In this case, Gq(x) — Gi(x) = exp[zi(a;— 1)]. In 
Eq. 0, <P k (R) = r[fc-l,zi(l-i?)]/r(fc - 1), where 
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F(n, x) is the incomplete gamma function. From Eq. (JTJ 
we get the size of the fc-core: 



M(fc) = p{l - T[k, *i(l - R)]/T(k)}, 



(11) 



where i? is the solution of Eq. 
gree distribution in the fc-core 



Q. The de- 
is P k (q>k) = 



pzl(l-R) q e~ Zl ^ R "^ /[M(k)q\]. Our numerical calcula- 
tions revealed that at p = 1, the highest fc-core increases 
almost linearly with z%, namely, kh ~ 0.78zi at z\ < 500. 
Furthermore, the mean degree z\(k) in the fc-core weakly 
depends on k: z\{k) « z\. 

Fig. |3 shows the dependence of the size of the fc-cores, 
M(k), on the concentration Q = 1 — p of the vertices 
removed at random. Note that counterintuitively, it is 
the highest fc-core — the central, most interconnected part 
of a network — that is destroyed primarily. The inset of 
Fig. [3 shows that with increasing damage Q, the mean 
degree Z\(k) decreases. The fc-cores disappear consecu- 
tively, starting from the highest core. The fc-core struc- 
ture of the undamaged ER graphs is displayed in Fig. 
Scale-free networks. — We consider uncorrelated net- 
works with a degree distribution P(q) oc (q + c)~ 7 . Let 
us start with the case of 7 > 3, where z 2 is finite. It 
turns out that the existence of fc-cores is determined by 
the complete form of the degree distribution including its 
low degree region. It was proved in Ref. that there is 
no fc ^ 3-core in a graph with the minimal degree q® = 1 , 
7 > 3, and c = 0. We find that the fc-cores emerge as 
c increases. The fc-core structure of scale-free graphs is 
represented in Fig. [3] The relative sizes of the giant fc- 
cores in the scale-free networks are smaller than in the 
ER graphs. As z 2 is finite, the k > 3-core percolation at 
7 > 3 is the hybrid phase transition. This is in contrast 
to the ordinary percolation in scale- free networks, where 
behavior is non-standard if 7 < 4 . 




FIG. 2: The size of the fc-core, M(fc), in the Erdos-Renyi 
random graph with the mean degree zi = 10 versus the con- 
centration Q of vertices removed at random. The highest core 
disappears at a very low concentration Q ~ 1.2% in contrast 
to the ordinary percolation threshold Q ~ 90%. The inset 
shows the mean degree 21 (fc) of vertices in the fc-core. 



The case 2 < 7 ^ 3 is realized in most important 
real-world networks. With 7 in this range, z 2 diverges if 
N — > 00. In the leading order in 1 — R -C 1, Eq. (J5J gives 
f k (R) (<j /fc) 7-2 (l - H)~ (3_7) . From Eq. © we nnd 
the order parameter R. Substituting this solution into 
Eq. (JIJ, in the leading order in 1 — R we find that the 
size of the fc-core decreases with increasing k: 



M(k)=p[q (l-R)/k 



7-1 



P 



2/(3- 



- 7 Wfc) (7 - 1)/(3 ~ 7) - 
(12) 

The divergence of fk{R) at R — > 1 means that the 
percolation threshold p c (k) tends to zero as N — > 00. 
The fc-core percolation transition in this limit is of infi- 
nite order similarly to the ordinary percolation [Icj . As 
kh{N — > 00) — ► 00, there is an infinite sequence of suc- 
cessively enclosed fc-cores. One has to remove at random 
almost all vertices in order to destroy any of these cores. 

Eq. ifTUj l allows us to find the degree distribution of 
fc-cores in scale-free networks. For 7 > 2 and fc 3> 1, 
Pk{q S> fc) ~ (7 — l)fc 7_1 g~ 7 . The mean degree z\{k) 
in the fc-core grows linearly with fc: zi(fc) ~ fczi/qo in 
contrast to the Erdos-Renyi graphs. 

Finite-size effect. — The finiteness of the scale-free net- 
works with 2 < 7 < 3 essentially determines their fc- 
core organization. We introduce a size dependent cutoff 
Qcut(N) of the degree distribution. Here q C ut(N) depends 
on details of a specific network. For example, for the con- 
figuration model without multiple connections, the de- 
pendence q cu t(N) ~ V^V" is usually used if 2 < 7 < 3. It 
is this function that must be substituted into Eqs. JT3J, 
((HJl . and below. A detailed analysis of Eq. (J5J 
shows that the cutoff dramatically changes the behavior 
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FIG. 3: The relative sizes of the fc-cores, M(k), panel (a), 
and fc-shells, iS'(fc), panel (b), in the Erdos-Renyi graphs with 
z\ — 10 and 20; scale-free networks with 7 = 2.5, 4 and 7, 
and an uncorrelated network with the degree distribution of 
the router- level Internet map (IR). The minimum degree in 
the scale-free networks is qo = 1. In the case 7 = 2.5, the 
maximum degree in the network is g cut = 2 000, and c = 2; 
for 7 = 4 and 7, c = 30 and 50, respectively. 
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of the function fk(R) near R=l. fk(R) has a maximum 
at i? max — 1— (3 — 7) _1 /(t -2 ) k/q cut and tends to zero at 
R — ► 1 instead of divergence. As a result, the fc-core per- 
colation again becomes to be the hybrid phase transition. 
The cutoff determines the highest fc-core: 



k h s p( 7 - 2) (3 - 7 ) (3 " 7)/(7 " 2) gcut(go/gc 



(13) 



The sizes of the fc-core at qo "C k <C fc^ are given by 
Eq. i|12|) . The relative size of the highest fc-core is 



M(fc fc )Sp[(3- 7 ) 



-(7-1)7(7-2). 



•l](9o/9e 



7 ~ 



Finally, the threshold of the fc-core percolation is 

Pc(k) = l//fe(i?max) = k/k h . 



(14) 



(15) 



If k—fkh, then p c (k) — > 1, i.e. even minor random dam- 
age destroys the highest fc^-core. By using exact Eqs. © 
and Q, we calculated numerically M(k) and S(k) for a 
scale-free network with 7 = 2.5, see Fig. [3] These curves 
agree with asymptotic expressions (fT^l and l(Tij). 

k-core organization of the router-level Internet. — We 
consider the router-level Internet which has lower degree- 
degree correlations than the Internet at the Autonomous 
Systems (AS) level. We substitute the empirical degree 
distribution of the router-level Internet as seen in skitter 
and iffinder measurements |l9l| into our exact equations 
and compare our results with the direct fc-core decompo- 
sition of this network. The calculated sizes of fc-cores and 
fc-shells are shown in Fig. |3| The calculated dependence 
S(k) [Fig.Efb), the IR curve] is surprisingly similar to the 
dependence obtained by the direct fc-core decomposition 
of, actually, a different network — the AS-level Internet — 
in Ref. [||. On the other hand, one can see in Fig.[3]that 
the highest fc-core with kh = 10 occupies about 2% of the 
network, while a direct fc-core decomposition of the same 
router-level Internet map in Ref. |]| revealed fc-cores up 
to kh = 32. This difference indicates the significance of 
degree-degree correlations, which we neglected. 

Discussion and conclusions. — It is important to indi- 
cate a quantity critically divergent at the fc-core's birth 
point. This is a mean size of a cluster of vertices of the fc- 
core with exactly fc connections inside of the fc-core. One 
may show that it diverges as —dM(k)/dp ~ (p — pc)" 1 / 2 
and that the size distribution of these clusters is a power 
law at the critical point. 

One should note that the fc-core (or bootstrap) perco- 
lation is not related to the recently introduced fc-clique 
percolation despite of the seemingly similar terms. 
The fc-clique percolation is due to the overlapping of fc- 
cliques — full subgraphs of fc vertices — by fc — 1 vertices. 
Therefore, the fc-clique percolation is impossible in sparse 
networks with few loops, e.g., in the configuration model 
and in classical random graphs, considered here. 

In summary, we have developed the theory of fc-core 
percolation in damaged uncorrelated networks. We have 



found that if the second moment of the degree distribu- 
tion of a network is finite, the fc-core transition has the 
hybrid nature. In contrast, in the networks with infinite 
Z2, instead of the hybrid transition, we have observed an 
infinite order transition, similarly to the ordinary perco- 
lation in this situation. All fc-cores in these networks are 
extremely robust against random damage. It indicates 
the remarkable robustness of the entire fc-core architec- 
tures of infinite networks with 7 < 3. Nonetheless, we 
have observed that the finite networks are less robust, 
and increasing failures successively destroy fc-cores start- 
ing from the highest one. Our results can be applied to 
numerous cooperative models on networks: a formation 
of highly connected communities in social networks, the 
spread of diseases, and many others. 
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